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^ The form factors of weak decays of the Be meson to orbitahy excited charmonium, 

^ , D, Bg and B mesons are calculated in the framework of the QCD-motivated rela- 

tivistic quark model based on the quasipotential approach. Relativistic effects are 
systematically taken into account. The form factor dependence on the momentum 
transfer is reliably determined in the whole kinematical range. The form factors are 
o ■ expressed trough the overlap integrals of the meson wave functions which are known 

J-., from the previous mass spectra calculations within the same model. On this basis 

D I semileptonic and nonleptonic B^ decay rates to orbitally excited heavy mesons are 

'"^ ' calculated. Predictions for the Be decays to the orbitally and radially excited 2P 

and 35 charmonium states are given which could be used for clarifying the nature 
^ ■ of the recently observed charmonium-like states above the open charm production 

G^ ! threshold. 
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PACS numbers: 13.20.He, 12.39.Ki 

I. INTRODUCTION 

The investigation of weak decays of mesons composed of a heavy quark and antiquark 
gives a very important insight in the heavy quark dynamics. The decay properties of the Be 
^ • meson are of special interest, since it is the only heavy meson consisting of two heavy quarks 

with different flavor. This difference of quark flavors forbids annihilation into gluons. As 
a result, the excited B^ meson states lying below the BD meson threshold undergo pionic 
or radiative transitions to the pseudoscalar ground state which is considerably more stable 
than corresponding charmonium or bottomonium states and decays only weakly. The CDF 
Collaboration reported the discovery of the Be ground state in pp collisions already more 
than ten years ago [l| . However, up till recently its mass was known with a very large error. 
Now it is measured with a good precision in the decay channel Be — )■ J/ipif. The measured 
value M*^^ = 6275.2 ± 2.9 ± 2.5 MeV [2| is in a very good agreement with the prediction of 
the relativistic quark model M^^""^ = 6270 MeV [3|]. More experimental data on masses and 
decays of the Be mesons are expected to come in near future from the Tevatron at Fermilab 
and the Large Hadron CoUider (LHC) at CERN. 

The characteristic feature of the Be meson is that both quarks forming it are heavy and 
thus their weak decays give comparable contributions to the total decay rate. Therefore it is 
necessary to consider both the b quark transitions b ^ c,u with the c quark being a spectator 
and c quark transitions c ^ s,d with the b quark being a spectator. The former transitions 



lead to weak decays to charmonium and D mesons while the latter lead to decays to Bg 
and B mesons. The estimates [^ of the B^ decay rates indicate that the c quark transitions 
give the dominant contribution (~ 70%) while the h quark transitions and weak annihilation 
contribute about 20% and 10%, respectively. However, from the experimental point of view 
the Be decays to charmonium are easier to identify. Indeed, CDF and DO observed the Be 
meson and measured its mass analyzing its semileptonic and nonleptonic decays Be — t- J /iplu 
and Be ^ J/ipn QMM- 

In this paper we extend our previous investigation of Be properties |3|, Ig, ITJ to study ex- 
clusive weak semileptonic and nonleptonic decay channels to orbitally excited heavy mesons. 
For the calculations we use the same effective methods (g, |7| developed in the framework 
of the relativistic quark model based on the quasipotential approach for the Be decays to 
ground and radially excited states of charmonium, D, Bg and B mesons. Here weak decays 
to orbital excitations of these mesons, governed both by the b and c quark decays, are consid- 
ered. The weak decay matrix elements are parametrized by invariant form factors which are 
then expressed through the overlap integrals of the meson wave functions. The systematic 
account for the relativistic effects, including wave function transformations from the rest 
to the moving frame and contributions from the intermediate negative-energy states, allows 
to reliably determine the momentum transfer dependence of the decay form factors in the 
whole accessible kinimatical range. The other important advantage of our approach is that 
for numerical calculations we use the relativistic wave functions, obtained in the meson mass 
spectra calculations, and not some ad hoc parametrizations which were widely used in some 
previous investigations. The calculated form factors are then substituted in the expressions 
for the differential decay rates. 

The important distinction between weak Be decays, associated with the h and c quark 
decays, consists in the significant difference of the accessible kinimatical ranges. In the Be 
decays to the charmonium and D mesons the kinimatical range is considerably broader (by 
about an order of magnitude) than for decays to Bs and B mesons. As a result, many 
weak decays which are kinematically allowed in the former case are forbidden in the latter 
one. The kinematical suppression of semileptonic Be — )■ Bs{B)lv decays should be more 
pronounced for the decays to excited states than for the ground ones. The nonleptonic 
Be decays to an orbitally excited heavy meson and an energetic light meson can then be 
considered on the basis of the factorization approximation. The obtained predictions for the 
decay rates are compared with previous calculations^ which are based on different relativistic 



quark models |8HlO|. three-point QCD sum rules yjj] and light-cone QCD sum rules [12 . 

We also consider here weak semileptonic and nonleptonic Be decays to the highly excited 
2P and 35 charmonium states. These states are of special interest since in last years a 
number of new charmonium-like states above the open charm production threshold have 
been observed [131]. They include several unexpectedly narrow states, X(3872), X(3940), 
y(3940), Z(3930), y(4260), Z(4430), which interpretation is controversial. Some of them 
could be candidates for excited charmonia. Therefore experimental observation of such 
states in Be decays could help to clarify their real nature. 



II. RELATIVISTIC QUARK MODEL 

In the quasipotential approach a meson is described as a bound quark- ant iquark state 
with a wave function satisfying the quasipotential equation of the Schrodinger type 

where the relativistic reduced mass is 

^^ - eTT^ " UP ' ^^^ 

and El , E2 are the center of mass energies on mass shell given by 

^' = 2M ' ^' = 2M • (^) 

Here M = Ei + E2 is the meson mass, mi^2 are the quark masses, and p is their relative 
momentum. In the center of mass system the relative momentum squared on mass shell 

reads 

2 . , [M^ - (mi + m2)^][M^ - (mi - m^)^] 

The kernel V{p, q; M) in Eq. ([1]) is the quasipotential operator of the quark-antiquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by [3|] 



V{p, q; M) = Ui{p)u2{-P)^{p, Q! M)ui{q)u2{-q), (5) 

with 

V(p, q; M) = ^asD,MlH + K%{k)TlT2;, + KL (k), 

where a^ is the QCD coupling constant, D^i, is the gluon propagator in the Coulomb gauge 

4-7T A-K ( k^k^\ 

D'%k) = -^, DH^) = -^ [S^' - ^j , D'^ = D^' = 0, (6) 

and k = p — q. Here 7^ and u{p) are the Dirac matrices and spinors 

u\p) 






where cr and x'^ cl'^g Pauli matrices and spinors and e{p) = y/p'^ + m^. The effective long- 
range vector vertex is given by 

^/.(k)=7M + ^^M-^^ (8) 



where k is the Pauh interaction constant characterizing the long-range anomalous chromo- 
magnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 

Vv{r) = {l-e){Ar + B), 

Vs{r) = e{Ar + B), (9) 

reproducing 

Konf(r) = Vs{r) + W(r) = Ar + B, (10) 

where e is the mixing coefficient. 

The expression for the quasipotential of the heavy quarkonia, expanded in v'^/c^ can be 
found in Ref. |3|. The quasipotential for the heavy quark interaction with a light antiquark 
without employing the nonrelativistic (v/c) expansion is given in Ref. p^]. All the param- 
eters of our model like quark masses, parameters of the linear confining potential A and 
B, mixing coefficient e and anomalous chromomagnetic quark moment k are fixed from the 
analysis of heavy quarkonium masses and radiative decays |3|. The quark masses mi, = 4.88 
GeV, rric = 1.55 GeV, m^ = 0.5 GeV, rriu^d = 0.33 GeV and the parameters of the linear 
potential A = 0.18 GeV^ and B = —0.30 GeV have values inherent for quark models. The 
value of the mixing coefficient of vector and scalar confining potentials e = —1 has been de- 
termined from the consideration of the heavy quark expansion for the semileptonic B ^^ D 
decays [l5| and charmonium radiative decays [3|. Finally, the universal Pauli interaction 
constant k = — 1 has been fixed from the analysis of the fine splitting of heavy quarko- 
nia ^Pj- states j3| and the heavy quark expansion for semileptonic decays of heavy mesons 



15l | and baryons [16|. Note that the long-range magnetic contribution to the potential in 
our model is proportional to (1 + k) and thus vanishes for the chosen value of k = — 1 in 
accordance with the fiux tube model. 



III. MATRIX ELEMENTS OF THE ELECTRO WEAK CURRENT FOR b ^ c,u 

AND c^ s,d TRANSITIONS 

In order to calculate the exclusive semileptonic decay rate of the B^ meson, it is necessary 
to determine the corresponding matrix element of the weak current between meson states. 
First we consider the weak B^ decays governed by the h quark decays. In the quasipotential 
approach, the matrix element of the weak current J^ = ^7^(1 — 75)6, associated with the 
6 — )■ g (g = c or m) transition, between a Be meson with mass Mb^ and momentum pb^ and 
a final P-wave meson F (F = Xcj, he or Dj ) with mass Mp and momentum pp takes the 
form [iTi] 

{F{pp)\J^\Be{pBj) = /^^^Fp,(p)r,(p,q)vl/^,p,^(q), (11) 

where r^(p, q) is the two-particle vertex function and "^mpm ^^^ the meson (M = Bc,F) 
wave functions projected onto the positive energy states of quarks and boosted to the moving 
reference frame with momentum Pm- 

The contributions to F come from Figs. [T] and [2J The contribution F*^^^ is the consequence 
of the projection onto the positive-energy states. Note that the form of the relativistic 
corrections emerging from the vertex function F^^^ explicitly depend on the Lorentz structure 
of the quark-antiquark interaction. In the leading order of the v'^/c^ expansion for Be and 
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FIG. 1: Lowest order vertex function F^^' contributing to the current matrix element (llip . 
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FIG. 2: Vertex function F'^^ taking the quark interaction into account. Dashed hues correspond to 
the effective potential V in ([5]). Bold lines denote the negative-energy part of the quark propagator. 



Xj and in the heavy quark limit rric — )■ oo for Dj only F*^^) contributes, while F*^^) contributes 
at the subleading order. The vertex functions look like 



F«(p, q) = Ugip.h^il - -f')ub{qb)i27f)'Sipr - qc), 



(12) 



and 






eb{k) + eb{pq) 



7?V(pc - qc 



+V(pc - qc 



eg(^') + ^qiQb) 



7?7im(1 - ll)}M(lb)uc{qc), 



(13) 



where the superscripts "(1)" and "(2)" correspond to Figs. [T] and [2], k = p^ — A; k' 
qfe + A; A = p^ - p^j 



A^~\p) 



e{p) - (m7° + 7°(7p)) 
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g.,c = 6,c(g)^±En«(m)g\ 





TABLE I: Mixing angles Lp for heavy-light mesons (in °). 




State 


D Ds B 


Bs 


IP 
2P 


35.5 34.5 35.0 
37.5 37.6 37.3 


36.0 
34.0 



and n*^*^ are three four-vectors given by 



i*r)J 



n 



^^^"(P)= 1^' '^^. + ir77^^^^ E = Jv' + M\ 



M' ' M{E + M) 
The wave function of a final P-wave F meson at rest is given by 

vl>^(p) ^ v]>^>',^,^^^(p) = yJM ^^(,,^,^^^(p), (14) 

where J and A4 are the total meson angular momentum and its projection, while S" = 0, 1 
is the total spin. ipp(^2s+ipj-j{p) is the radial part of the wave function, which has been 



determined by the numerical solution of Eq. ([T]) in |3|, [ij]. The spin-angular momentum 
part ys''^ has the following form 

ys^= E(l-^-^i-^2, Sai+a2\JM){lau Ict2\S ai+a2)Y{^-'''-'''xMi)X2{cT2). (15) 

Here {jirrii, J2^^2\J M) are the Clebsch-Gordan coefficients, Y[^ are spherical harmonics, 
and x(c) (where a = ±1/2) are spin wave functions, 

X(1/2)=Q, x(-l/2)=(°). 

The heavy-light meson states (such as Di, D[ etc.) with J = L = 1 are mixtures of 
spin-triplet F(^Pi) and spin-singlet F(^Pi) states: 

"^Fi = '^F(^Pi) cos (p + \I/f(3Pi) sin (f, 

-qtp^ = -'^p^^ip^)sm(p + '^F(3Pi)Cos(p, (16) 

where (p is the mixing angle and the primed state has the heavier mass. Such mixing occurs 
due to the nondiagonal spin-orbit and tensor terms in the Qq quasipotential. The physical 
states are obtained by diagonalizing the corresponding mixing terms. The values of the 
mixing angle (p were determined in the heavy-light meson mass spectra calculations [l4'] and 
are given in Table [11 

It is important to note that the wave functions entering the weak current matrix element 
fITT]) are not in the rest frame in general. For example, in the B^ meson rest frame (p^^ = 0), 
the final meson is moving with the recoil momentum A. The wave function of the moving 
meson \1/fa is connected with the wave function in the rest frame "^fo = ^f by the 
transformation |17| 

^fa(p) = Dl/\RYjDl/\RYj^Fo{p), (17) 

where R^ is the Wigner rotation. La is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D^^'^(R) in spinor representation is given by 

{l ^^Dli\RYJ = S-\p,,^)S{A)Sip), (18) 



where 



Q. X /e(p)+m / ctp 



is the usual Lorentz transformation matrix of the four-spinor. 

The expressions for the matrix elements of the B^ decays to the P-wave Bgj and Bj 
mesons, governed by the c quark decays, can be obtained from the above expressions by the 
interchange of the b and c quarks and for the final active quarks q = s,d. 

IV. FORM FACTORS OF THE SEMILEPTONIC B^ DECAYS TO THE 
ORBITALLY EXCITED HEAVY MESONS 

The matrix elements of the weak current J^ = &7/i(l — 75)5 or 07^ (1 — 75)5 for Be decays 
to orbitally excited scalar light mesons [S) can be parametrized by two invariant form factors 

{S{pF)\qYb\B,{pBj) = 0, 

{SiPF)\qi''l,b\B,{pBj) = U{q') {p'b. + p'f) + /-(?') {pk " Pf) , (19) 

where q = pb^ — Pf, Ms is the scalar meson mass. 

The matrix elements of the weak current for B^ decays to axial vector mesons (AV) can 
be expressed in terms of four invariant form factors 

{A{pF)\qrb\B,{pBj) = {Mb^ + MA)hy,{q')e*^ + [/^^.(g'K. + hy.iqV^f-^, (20) 

Mb, 

{A{pF)\qri,b\B^{pBj) = ,!'^^it] e''''''<PB.,PF., (21) 

Mbc + ^vIa 

where Ma and e'' are the mass and polarization vector of the axial vector meson. 

The matrix elements of the weak current for Be decays to tensor mesons (T) can be 
decomposed in four Lorentz-invariant structures 

(T(Pf)|«i''75()|B=(pbJ) = (Afe, + Afr)(A,(«')e-'-^ 

+ [tM{e)pk + ^a^VfYJ^^, (23) 

where Mt and e^'^ are the mass and polarization tensor of the tensor meson. 

We previously studied the form factors parametrizing the matrix elements of vector and 
axial vector charged and neutral weak currents for Be — ?■ r]c{J/4'), Be — )■ -D*-*-* [6], Be — )■ 
B^*\B^*^) [7| and Be — > -D^*-* transitions in the framework of our model. Now we apply 
the same approach, described in detail in Refs. [a, [7|, [l8|, for the calculation of the form 
factors for Be decays to the orbitally excited heavy mesons. Namely, we calculate exactly 
the contribution of the leading vertex function T^^^ flT2l) to the transition matrix element of 
the weak current flTTl) using the (5-function. For the evaluation of the subleading contribution 
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r^^) for the Be — )■ Xj{hc) and Be — > Dj transitions, governed by 6 — )■ c, m transitions, we 
use expansions in inverse powers of the heavy 6-quark mass from the initial Be meson and 
large recoil energy of the final heavy meson. Note that the latter contributions turn out to 
be rather small numerically. Therefore we obtain reliable expressions for the form factors 
in the whole accessible kinimatical range. It is important to emphasize that doing these 
calculations we consistently take into account all relativistic corrections including boosts 
of the meson wave functions from the rest reference frame to the moving ones, given by 
Eq. (fTTl) . The obtained expressions for the decay form factors are given in Appendix (to 
simplify these expressions the long-range anomalous chromomagnetic quark moment was 
explicitly set as k = — 1). In the limits of infinitely heavy quark mass and large energy of 
the final meson, the form factors in our model satisfy all heavy quark symmetry relations 



m 20 



As a result, we get the following expressions for the Be decay form factors: 

(a) Be ^ S transitions {S = Xco, Dq) 

Uq') = /±\q')+ef^^'\q') + {1 - s)fl^'\q% (24) 

(b) Be -^ AV transition {AV = Xci, /^i(^Pi)) 

hvAQ') = h^v^q') + ehf\q') + {I - e)hl^\q') , {i = 1,2,3), 

hA{q') = h^hq') + ehf\q') + {I ~ e)h^/\q\ (25) 

(c) Be -^ AV transition^ {AV = he,DiCPi)) 

9vAq') = 9^vX<f) + e9v'\<l') + {l-e)gl!^\q'), (^ = 1'2,3), 

9A{q') = 9^^{q') + e9T\q') + (1 " e)gT{q'). (26) 

(d) Be^T transition (T = Xc2, ^2) 

tv{q') = t^v\q') + e4^'\q') + {I - e)tl^'\q\ 

tA^q') = t^Z{q') + etf\q^) + {l-e)t'2^\q% {^ = 1,2,3), (27) 

where fl^^ ff'^^^^ h^^'^ h^'^^"^^ h^^^ h^'^^^^ n^^^ n^'^^^^ n^^^ o^'^^^^ /^^ /^'^^^^ f^^^ and 
wnere j± ■, j± , iiy. , ay. , n^^ , n^^ , gy. , gy. , g^^ , g^^ , ly , ly , z;^^ , ana 

t^'. are given in Appendix. The superscripts "(1)" and "(2)" correspond to Figs. [1] and [21 
S and V correspond to the scalar and vector confining potentials of the gg-interaction. The 
mixing parameter of scalar and vector confining potentials e is fixed to be —1 in our model. 
In the case of Be decays to P-wave Bg and B mesons, governed by the c ^ s,d transitions, 
the accessible kinematical range is significantly smaller (by almost a factor of 4) than the 
one for the decays to the S'-wave Bg and B mesons. Our previous investigation [7| of the 
latter decays had shown that intermediate negative-energy states, leading to the subleading 
term T^'^\ give almost negligible contributions to decay form factors (see Fig. 3 of Ref. |7|). 



^ The corresponding decay matrix elements are defined by Eqs. (|20|) and (|2T|) with the replacement of form 
factors hi[q^) by gi{q^)- 



Therefore such contributions can be safely neglected in the present analysis. Thus, for 
calculations of the form factors of the B^. — )■ Bgj and B^ — )■ Bj weak transitions we use 
the leading order expressions /j , h] , gl and t} , given in Appendix, where the b and c 
quarks are interchanged. 

For numerical calculations of the form factors we use the quasipotential wave functions 
of the Be meson and orbitally excited charmonium and D, Bs, B mesons obtained in their 
mass spectra calculations |3|, [iJ] . Our results for the masses of these mesons are in good 



agreement with available experimental data ^]. Therefore we use the experimental values 
for the masses of well-established states and our model predictions for all other masses in 
the numerical calculations. 

In Fig. [3] we plot form factors of the Be weak transitions to the IP {XcJ, he) and 2P 
{x'cjj h'c) -wave charmonium states as an example. The remaining plots for the B^ weak 
form factors to the P-wave Dj, Bgj and Bj mesons have an analogous behaviour and are 
not shown here. 



V. SEMILEPTONIC B^ DECAYS TO ORBITALLY EXCITED HEAVY MESONS 

The differential decay rate for the B^ meson decay to P-wave heavy mesons reads [Sj 



dT{Be-^F{S,AV,T)W) Gl .^^ ,2>^^'\q' - 



H, 



m 



2\2 



™*(-S)-|?-'-.' 



rfg2 (27r)3' ^^' 24M|g2 

(28) 
where Gp is the Fermi constant, Vij are the Cabbibo-Kobayashi-Maskawa (CKM) matrix 
elements, A = X{Ml^,Mj,, q^) = M%^ + M^ + q'^ - 2(M|^M|, + Mj,q'^ + M^j'^), mi is the 
lepton mass and 

HH^ = H+hI + H_h1_ + HoHI (29) 

Helicity components of the hadronic tensor are expressed through the invariant form factors, 
(a) Be — !■ S(^Po) transition 



H± = 

Ai/2 

3 



Ho = ^^U{q^) 



Ht = -^ml - MDUiq') + q'nq% (30) 

(b) Be -^ AV(^Pi) transition 
H^ = {Mb^ + MAv)hvAQ^) ± 7^-^ri^— ^A, 

Mb,^ + J-vJ-AV 

Ho = WTT^ {{Mb^ + MavKMI - Mly - q')hyM') + l^[hv.{q') + hvM')]] , 
2MAv\/q^ { ^Mb, J 

Ht = \{MB^ + MAy)hy,{q')+ ^; ^"^ [hv,{q') + hy,{q')] 

2MAvVq [ 2Mb, 

-^l^[hvAq')'hyM')]\ (31) 
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FIG. 3: Form factors of the B^ decays to the IP- and 2P -wave charmonium states. 
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(c) Be -^ AV'i^Pi) transition 

Hi are obtained from expressions (!3T!) by replacement of form factors hi{q^) by gi{q^)- 

(d) B^ -^ T(3P2) transition 



H 



Ai/2 



(M.. + M.)t.,(,^)±^^-^t. 



"' ■[tA.(g^)-t^3 (?')])• (32) 



2Mf 



Here the subscripts ±,0,t denote transverse, longitudinal and time helicity components, 
respectively. 

Now we substitute the weak decay form factors calculated in the previous section in the 
above expressions for decay rates. The resulting differential distributions for B^ decays to 
the IP (xj, he) and 2P ix'ji h'J charmonium states are plotted in Fig. HI The difference of 
the plot shapes for the corresponding IP and 2P charmonium states is the consequence of 
their different nodal structure. We calculate the total rates of the semileptonic Be decays 
to the P-wave heavy mesons by integrating the corresponding differential decay rates over 
g^. For calculations we use the following values of the CKM matrix elements: \Vcb\ = 0.041, 
\Vub\ = 0.0038, \Vcs\ = 0.974, \Vcd\ = 0.223. Our predictions for the rates of the semileptonic 
Be decays to the P-wave charmonium states are compared with the previous calculations 
[8|-[l2] in Table HTl The authors of Refs. |8l-ll0| use different types of relativistic quark models. 



Calculations in Ref. pj| are based on the three-point QCD sum rules, while Ref. [1^ employs 
the light-cone QCD sum rules. We find that significantly different theoretical approaches 
give values for the B^ — )■ Xj{h^lv decay rates consistent in the order of magnitude, while 
for the Be decays to first radial excitations of the P-wave charmonium {Be — ?■ 'x'j{h'^)lv) 
our results are almost an order of magnitude lower than the predictions of the light-cone 
QCD sum rules [12], which are the only available ones at present. The latter decays can 
play an important role in studying charmonium states above the open charm production 
threshold. Their observation at Tevatron and LHC can help to clarify the nature of the 
new charmonium-like states. Our results for the rates of the CKM suppressed semileptonic 
Be decays to the P-wave D mesons, governed by the weak h ^ u transitions, are given in 
Table 

In Fig. O we plot predicted differential semileptonic decay rates of the Be to P-wave Bg 
meson states, governed by the c — )■ s weak transitions. The allowed kinimatical range for 
these transitions is rather narrow. Therefore semileptonic decays involving the r lepton are 
forbidden. From these plots we see that even the account of the rather small muon mass 
significantly modifies the differential decay rates. The corresponding plots for Be — ?■ Bjlu 
decays have similar shape and are not shown here. The predicted values for the rates of the 
semileptonic Be decays to the P-wave Bg and B mesons are given in Tables [IV] and El Note 
that, notwithstanding the fact that these decay rates have significantly larger values of the 
CKM matrix elements than the rates of Be decays to charmonium, they have the same order 
of magnitude. This is the result of the above mentioned strong phase space suppression of 
the Be — !■ Bsjlv decays. 
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FIG. 4: Predictions for the differential decay rates of tlie Be semileptonic decays to the IP- and 
2P-wave charnionium states. 
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TABLE II: Comparison of our predictions for the rates of the semileptonic Be decays to the P-wave 
charmonium states with previous calculations (in 10^^^ GeV). 



Decay 


our 


[8] 


[9] 


[10] 


[11] 


[12] 


Be -^ Xcoeu 


1.27 


2.52 


1.55 


1.69 


2.60 ±0.73 




Be -^ XcQTU 


0.11 


0.26 


0.19 


0.25 


0.70 ± 0.23 




Be -^ Xciev 


1.18 


1.40 


0.94 


2.21 


2.09 ± 0.60 




Be -^ XclTf 


0.13 


0.17 


0.10 


0.35 


0.21 ±0.06 




Be -^ Xc2eiy 


2.27 


2.92 


1.89 


2.73 






Be -^ Xc2Tl' 


0.13 


0.20 


0.13 


0.42 






Be — )■ heCU 


1.38 


4.42 


2.40 


2.51 


2.03 ±0.57 


4.2 ±2.1 


Be — )■ hcTU 


0.11 


0.38 


0.21 


0.36 


0.20 ±0.05 


0.53 ±0.26 


Be -^ X'eo^^'^ 


0.19 










10 ±6 


Be -^ x'eQT^ 


0.0089 










0.39 ±0.20 


Be -^ Xcie^ 


0.12 










8.6 ±4.8 


Be -^ X'el^l^ 


0.0056 










0.31 ±0.18 


Be -^ Xc2e^ 


0.048 












Be -^ Xc2^'^ 


0.0019 












Be -)• h'eev 


0.031 










0.76 ± 0.33 


Be -)• h'eTV 


0.0016 










0.028 ± 0.014 



TABLE III: Predictions for the rates of the semileptonic Be decays to the P-wave D mesons (in 
10-15 GeV). 



Decay 



Decay 



Be -> D^ev 
Be — ^ Diev 
Be -^ D[ev 



Br 



Doev 



0.016 
0.016 

0.027 
0.052 



Br 



DoTU 



Be — ^ DiTV 
Be -^ D[tU 



Br 



DoTV 



0.0067 
0.0056 
0.016 
0.019 



VI. NONLEPTONIC DECAYS 



In the standard model nonleptonic Be decays are described by the effective Hamiltonian, 
obtained by integrating out the heavy VT-boson and top quark, 
(a) For the case of the 6 — )• c, m transitions, one gets 



-f^eff — ~~7K^cb 



ci(/i)0f±C2(/x)0: 



eb 



V2 



K 



ub 



Ci(/i)Of ±C2(/i)0^ 



(b) For the case of the c — )■ s, d transitions, we have 

H,s = ^Ves [Cl(/i)Or ± C2(/i)0r] ± ^Ved [ci(/i)OJ'^ ± C2(/i)0^'^ 



± 



+ 



(33) 



(34) 



The Wilson coefficients Ci^2(Ai) are evaluated perturbatively at the W scale and then are 
evolved down to the renormalization scale /j, ~ nib by the renormalization-group equations. 



14 




0.10 

g2 (GeV^) 



0.15 



0.20 



B,^B,,lv, 




q^ (GeV^ 



0.000 
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FIG. 5: Predictions for the differential decay rates of the Be semileptonic decays to the P-wave Bg 
meson states. 



TABLE IV: Predictions for the rates of the semileptonic Be decays to the P-wave Bg mesons (in 
10-15 GeV). 



Decay 



Decay 



Be -> Bsoeu 
Be -^ Bsieu 
Be -^ B'^^eu 



Be 



Bs2ev 



0.96 
0.029 
0.065 
0.066 



Be — )• BsQ^xy 

Be -^ Bsi^JLV 
Be -^ B'^^fll^ 



Br 



Bs2fJ'l' 



0.82 
0.026 
0.044 
0.031 



TABLE V: Predictions for the rates of the semileptonic Be decays to the P-wave B mesons (in 
10-15 GeV). 



Decay 



r 



Decay 



r 



Be — )• Boeu 
Be — ^ Bieu 
Be -^ B[ev 



Br 



B2ev 



0.089 
0.0048 
0.010 
0.012 



Be -^ BoflU 
Be — )• Bi^V 
Be -^ B'lflU 



Br 



B2IJ.I' 



0.082 
0.0043 
0.0082 
0.0067 
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The ellipsis denote the penguin operators, the Wilson coefficients of which are numerically 
much smaller than Ci^2- The local four-quark operators Oi and O2 are given by 

of = [idu)v-A + isc)v-A]iqb)v-A, 

of = {qu)v-A{dh)v-A + {qc)v-A{sh)v-A, q = {u,c), (35) 

and 

Oi'' = {du)v-A{cq)v~A, 

02^ = {cu)v-A{dq)v-A, q=is,d), (36) 

where the rotated antiquark fields are 

d = Vudd + KsS, s = Vcdd + VcsS, (37) 

and for the hadronic current the following notation is used 

{qq')v-A = qit^{i - i5)q' = J^ ■ 

The factorization approach, which is extensively used for the calculation of two-body 
nonleptonic decays, such as B^ — )■ FM, assumes that the nonleptonic decay amplitude 
reduces to the product of a meson transition matrix element and a decay constant |22| . 
This assumption in general cannot be exact. However, it is expected that factorization 
can hold for energetic decays, where the final F meson is heavy and the M meson is light 
A justification of this assumption is usually based on the issue of color transparency 



23| 

24j |. In these decays the final hadrons are produced in the form of almost point-like color- 
singlet objects with a large relative momentum. And thus the hadronization of the decay 
products occurs after they are too far separated for strongly interacting with each other. 
That provides the possibility to avoid the final state interaction. A more general treatment 
of factorization is given in Refs. |25|, |26| . 



Here we first analyze the B^ nonleptonic decays to the P-wave charmonium and the 
light TT"*", p+ or i^(*)+ mesons, governed by the weak 6 — )■ c, m transitions. The corresponding 
diagram is shown in Fig. EJ^a), where qi = d, s and g2 = u. Then the decay amplitude can 
be approximated by the product of one-particle matrix elements 

{F'M+\H,s\Bt) = ^V,,V,,,MF\{bc)v-A\Bc){M\{q,q2)v-A\0), (38) 



where 



oi = ci(/i) + — C2(/i) (39) 



'c 

and Nc is the number of colors. 

Next we consider nonleptonic decays of the Be meson to the P-wave Bg or B mesons 
and the final light M+ meson, governed by the weak c ^ s,d transitions. Only the pion 
is kinematically allowed. The corresponding diagram is shown in Fig. Mjo). Then in the 
factorization approximation the decay amplitude can be expressed through the product of 
one-particle matrix elements 

(P°M+|/feE|5+) = ^V,,Vg,,MF\{cq)v^A\B,){M\{qiq2)v^A\0). (40) 
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qi 



92 



M+ 




A•/^ 



5^ 



F' = Xcj, K 



i?^ 



F' = B,j,BJ 



(a) (b) 

FIG. 6: Quark diagram for the nonleptonic B^ — )• F^M^ decay. 

The matrix element of the weak current J^ between vacuum and a final pseudoscalar 
(P) or vector (V^) meson is parametrized by the decay constants fpy 



{P\qiY'y5q2\0) = ifpp^p, (y\q^lM^) = ^i^Myh 



(41) 



The pseu doscalar fp and vector fv decay constants were calculated within our model in 



Ref. [271]. It was shown that the complete account of relativistic effects is necessary to 



get agreement with experiment for decay constants especially of light mesons. We use the 
following values of the decay constants: f-,, = 0.131 GeV, fp = 0.208 GeV, fx = 0.160 GeV 
and Jk* = 0.214 GeV. The relevant CKM matrix elements are |Kd| = 0.975, |Ks| = 0.222. 
The matrix elements of the weak current between the B^ meson and the final heavy meson 
F entering the factorized nonleptonic decay amplitude (138!) are parametrized by the set of 
decay form factors defined in Eqs. fl2^ - (p71) . Using the form factors obtained in Sec. llVt we 



get predictions for the nonleptonic B^ — )■ yj(/icl°-^^ decay rates and give them in Table 



in comparison with other calculations |8l-llOL |28| . which are available for the decays to the 
IP charmonium states only. Predictions for the energetic nonleptonic decays to the 2P 
charmonium states are made for the first time and their measurement could be important 
for the identification of these states. Our results for the nonleptonic B'^ — )■ BsjM^ and 
5+ — )■ BjM^ decay rates are presented in Table IVIIi Note that in the latter case only 
decays involving pions are kinematically allowed. 



VII. CONCLUSIONS 



We calculated the form factors of the weak Be decays to orbitally excited heavy mesons, 
governed both by the b ^ c,u and c — )■ s, rf transitions, in the framework of the QCD- 
motivated relativistic quark model based on the quasipotential approach. The momentum 
dependence of the weak decay form factors was reliably determined in the whole accessible 
kinematical range. This is particularly important for Be decays to the P-wave charmonium 
and D mesons since they have a rather broad kinematically allowed range (g^ax ~ 6 — 
15 GeV^). All essential relativistic effects were taken into account including transformations 
of the meson wave functions from the rest to the moving reference frame and contributions 
from the intermediate negative-energy states. The resulting form factors are expressed 
through the overlap integrals of the meson wave functions. These wave functions were 
obtained previously in the meson mass spectra calculations and are used in the present 
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TABLE VI: The 


rates of the nonleptonic Be decays 


to the P-wave 


charmonium and hght 


mesons 


(in 10-15 GeV). 












Decay 


our 


[8] 


[9] 


[10] 


[28] 


B+ ^ XcO^+ 


0.23af 


0.622a? 


0.28a? 


0.317a? 


11a? 


Bt -^ XcOP^ 


0.64af 


1.47a? 


0.73a? 


0.806a? 


37a? 


B+ ^ XcoK+ 


0.018a? 


0.0472a? 


0.022a? 


0.00235a? 




B+ ^ XcoK*+ 


0.045a? 


0.0787a? 


0.041a? 


0.00443a? 




B+ ^ xcivr+ 


0.22af 


0.0768a? 


0.0015a? 


0.0815a? 


0.10a? 


i?c+ ^ XC1/0+ 


0.16af 


0.326a? 


0.11a? 


0.331a? 


5.2a? 


S+ ^ Xcli^+ 


0.016a? 


0.0057a? 


0.00012a? 


0.0058a? 




5+ ^ Xcli^*+ 


0.010a? 


0.0201a? 


0.0080a? 


0.00205a? 




5c+ ^ Xc2^+ 


0.41a? 


0.518a? 


0.24a? 


0.277a? 


8.9a? 


Bt ^ Xc2P+ 


1.18a? 


1.33a? 


0.71a? 


0.579a? 


36a? 


B+ ^ Xc2i^+ 


0.031a? 


0.0384a? 


0.018a? 


0.00199a? 




5+ ^ Xc2i^*+ 


0.082a? 


0.0732a? 


0.041a? 


0.00348a? 




5+ -^ hcTr+ 


0.51a? 


1.24a? 


0.58a? 


0.569a? 


18a? 


B+ -^ hcp+ 


1.11a? 


2.78a? 


1.41a? 


1.40a? 


60a? 


B+ -^ hcK+ 


0.039a? 


0.0939a? 


0.045a? 


0.0043a? 




B+ -^ hcK*+ 


0.077a? 


0.146a? 


0.078a? 


0.0076a? 




Bt ^ x:.o^+ 


0.023a? 










Bt ^ ^coP-^ 


0.080a? 










Bt ^ x'cO^+ 


0.0019a? 










Bt ^ x^o^*+ 


0.0055a? 










i?c+ ^ X^ivr+ 


O.OUa? 










i?c+ ^ XciP^ 


0.016a? 










i?+ ^ Xcl^+ 


0.0095a? 










5+ ^ Xcli^*+ 


0.0011a? 










Bt ^ X'c2^+ 


8.5xl0-^a? 










Sc+ ^ X^2P+ 


0.0022a? 










S+ ^ X'c2^+ 


8.5xl0-'^a? 










5+ ^ X^2^*+ 


0.00015a? 










S+ ^ /l>+ 


l.OxlO-^a? 










Bt -^ Kp+ 


0.0051a? 










Bt -^ KK+ 


3.4xl0-6a? 










Bt -^ hcK*+ 


0.00035a? 











numerical evaluations. The influence of mixing effects on the P-wave heavy-light meson 
wave functions due to the non-diagonal spin-orbit and tensor terms in the Qq quasipotential 
was explicitly considered. The reliable determination of the q^ dependence of the from 
factors in the whole kinematical range is an important achievement, since in many previous 
calculations form factors were determined only at the single point of either zero (g^ = Q'^ax) 
or maximum (g^ = 0) recoil of the final meson, and then different ad hoc extrapolations 
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TABLE VII: The rates of the nonleptonic B^ decays to the P-wave Bs or B mesons and vr meson 
(in 10-15 GeV). 

Decay F Decay T 

B+ -^ BsoTT+ 5.82af B+ -^ B^Tr+ 0.46a? 

B+ -^ BsiTT+ O.SOaj B+ -^ B^Tr+ OMlaj 

B+ -^ B'^^7r+ 0.31a? B+ -^ B[^Tr+ 0.061a? 

B+ -^ Bs2TT+ 0.26a? B+ -^ B^tt+ 0.047a? 



were employed. 

The obtained weak form factors were used for the calculation of the semileptonic and 
nonleptonic Be decays to corresponding orbitally excited heavy mesons. For the nonleptonic 
decays the factorization approximation was used. The calculated branching fractions are 
summarized in Table IVIIII In this table we give our predictions not only for Be decays 
to the first IP- wave charmonium states {xj,hc), but also for their radial excitations (2P- 
wave charmonium Xj,/i^). For completeness, we also present there our predictions for the 
semileptonic P^decays to the 35" charmonium states {ip",f]c) which were not given in our 
previous study |6[. These decays to highly (both radially and orbitally) excited charmonium 
are of special interest, since their observation could help to reveal the nature of the newly 
observed charmonium-like states above the open charm production threshold. 

Summing the corresponding branching fractions in Table IVIIII we find that the semilep- 
tonic^ and the considered energetic nonleptonic decays to the IP charmonium states con- 
tribute about 0.88% and 0.44% of the total rate, respectively. The corresponding decays to 
the 2P charmonium states are significantly suppressed (by an order of magnitude) mainly 
due to the presence of the node in the 2P wave function and give about 0.057% and 0.013% 
of the total rate. The same pattern was previously observed in B^ decays to the S-wave 
charmonia |6| , where the rates of decays to the 25* states were also suppressed by an order of 
magnitude compared to decays to the IS" states. For decays to higher charmonium excita- 
tions such suppression should be even more pronounced. The CKM suppressed semileptonic 
decays to the Dj mesons contribute about 0.019%. Thus the total contribution of the 
considered Be decays, governed by the b ^ c,u weak transitions, is about 1.41%. 

The Be semileptonic decays to orbitally excited Bgj and Bj mesons, governed by the 
c — )■ s, d weak transitions, turn out to have smaller branching fractions than Be decays to 
orbitally excited charmonium, notwithstanding the significantly larger values of the CKM 
matrix elements, due to the substantial phase space suppression. Such decays involving 
the T are kinematically forbidden, while the muon mass starts to play an important role, 
reducing branching fractions by more than 10%. In total, such semileptonic decays give 
about 0.045% of the Be decay rate. On the other hand, there is no kinematical suppression 
in the corresponding nonleptonic decays and they contribute about 0.69%. Thus, the total 
contribution of the considered Be decays, governed by the c — )■ s, d weak transitions, is about 
0.73%. 

The semileptonic and nonleptonic Be decays to excited heavy mesons can be investigated 



^ We also take into account semileptonic decays involving the muon, which rates are almost equal to the 
ones with the electron. 



TABLE VIII: Branching fractions (in %) of exclusive 
the Be hfetime tb^ = 0.46 ps and ai = 1.14 for the b — 
transitions. 
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Be decays calculated for the fixed values of 
■ c transitions and oi = 1.20 for the c ^ s,d 



Decay 


Br 


Decay 


Br 


Decay 


Br 


Be -^ Xcoei^ 


0.087 


Bt -^ Xcovr+ 


0.021 


Bt ^ X'eO^"- 


0.0020 


Be -^ XcoTiy 


0.0075 


Bt -^ XcoP+ 


0.058 


Bt ^ X'eOP^ 


0.0071 


Be -^ Xcieu 


0.082 


Bt ^ Xe^K+ 


0.0016 


Bt ^ x'eoK+ 


0.00017 


Be -^ XciTiy 


0.0092 


B+ ^ Xe^K*+ 


0.0040 


Bt ^ x'eoK*+ 


0.00049 


Be -^ Xc2ey 


0.16 


Bt -^ Xcivr+ 


0.020 


Bt ^ X>^ 


0.0010 


Be -^ Xc2Ty 


0.0093 


Bt -^ XelP+ 


0.015 


Bt -^ XelP+ 


0.0014 


Be — > hetU 


0.096 


Bt ^ XelK+ 


0.0015 


Bt ^ x'eiK^ 


0.000086 


Be — )■ heTU 


0.0077 


Bt ^ XelK*+ 


0.0010 


Bt ^ x'eiK*^ 


0.00010 


Be -^ x'eO^V 


0.014 


Bt -^ Xc2vr+ 


0.038 


Bt ^ X>^ 


7.7x10^8 


Be -^ x'eOT'^ 


0.00063 


Bt ^ Xe2P+ 


0.11 


Bt ^ X'e2P-^ 


0.00020 


Be -^ X'el^l^ 


0.0085 


Bt ^ Xe2K+ 


0.0028 


Bt ^ X'e2K^ 


7.8x10"^ 


Be -^ X'el^y 


0.00039 


Bt ^ Xe2K*-^ 


0.0074 


Bt ^ x'e2K*+ 


0.000014 


Be -^ Xc2eJ^ 


0.0033 


Bt -^ heTT+ 


0.046 


Bt -^ h'e7r+ 


9.4x10-^ 


Be -^ X'e2^y 


0.00013 


Bt -^ hep+ 


0.10 


Bt -^ Kp+ 


0.00046 


Be — > h'eev 


0.0021 


Bt -^ heK+ 


0.0035 


Bt -^ KK+ 


3.1x10"^ 


Be -^ h'eTU 


0.00011 


Bt -^ heK*+ 


0.0070 


Bt -^ h'eK*+ 


0.000032 


Be -^ Doeu 


0.0011 


Be -> Bsoei^ 


0.0066 


Be — )■ BqEU 


0.0061 


Be — ^ DoTl/ 


0.00046 


Be -^ BsQp.V 


0.0057 


Be — > Bq^V 


0.0056 


Be -^ Diev 


0.0011 


Be -^ Bsiev 


0.0020 


Be -^ Bieu 


0.00033 


Be -^ DiTU 


0.00039 


Be — ^ Bslj-liy 


0.0018 


Be -^ Bi^U 


0.00030 


Be -^ D[ev 


0.0019 


Be -^ B'^^eu 


0.0045 


Be -^ B[eu 


0.00072 


Be -^ D[tV 


0.0011 


Be — ^ B'g^^u 


0.0031 


Be -^ B[flU 


0.00057 


Be -^ D2eu 


0.0036 


Be -^ Bs2ey 


0.0046 


Be — > B2eu 


0.00084 


Be — > D2TU 


0.0013 


Be — > Bs2P'I^ 


0.0022 


Be -^ B2PV 


0.00047 


Be -^ r]'^ev 


0.00055 


Bt -^ BsQ^+ 


0.55 


Bt -^ BoTT+ 


0.043 


Be -^ i!,TU 


5.0x10^^^ 


Bt ^ Bsi^+ 


0.028 


Bt -^ Bi7r+ 


0.0039 


Be -^ ip"eu 


0.00057 


Bt ^ B',^7i+ 


0.029 


Bt ^ B'o7r+ 


0.0058 


Be -^ '4)"tU 


3.6x10-6 


Bt -^ Bs2TT^ 


0.024 


Bt -^ B2TT+ 


0.0044 



at Tevatron and LHC, especially in the LHCb experiment, where the Be mesons are expected 
to be copiously produced. 
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Appendix: Form factors of weak Be decays to orbitally excited heavy mesons 

a) Be -)■ S{^Po) transition {S = Xco, -Dq) 

2mr 



E, 



X 



+ 



Mb J (27r 
(pA) 



:i's{p + 



jip + A) +mg 



Es + Ms J\ 2e,(p + A) \ 2eb{p) 
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2 p 
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I 1 



±{Mb^TEs)\\ 
1 ^ 



P^ 



2 



eg(p + A) + mg] \ei,{p) + m^] / 
A2 



X 1 



P 



[eg(p + A) + mg][e6(p) + mj 

>^Bc(P 



+ P 
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1 

+ 



eg(p + A)+mq efe(p)+mb 
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X It 
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2mr 



1 



-p\ 



-^Ap + 



+ 
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Mb^ J {27if ^'V^^ Es + Ms )\ 2e,(A) 
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X- 



p 
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' 3eg{A)[eg{A) + mg 
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2mr 



-A 
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1 [^ ^ Mb^ t Es 



2mLh 



eg (A) +mg 
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Es + Ms J V Es + Ms 

(b) 5c ^ AV(3Pi) transition {AV = Xci,Di{^Pi)) 



-A 



V'Bc(p) 
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where the subscript q = c,u refers to the final active quark, the superscripts "(1)" and "(2)" 
correspond to Figs. [1] and [21 S and V correspond to the scalar and vector potentials of the 
gg-interaction, and 



A= lAI 



\ 



(Ml + Ml- g2)2 



4M|^ 



-Ml, {F = S,AV,AV\T) 



Ef 



M2+A2, eQ{p + aA) = Jml + {p + aAy {Q = b,c,s,u,d). 
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